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ABSTRACT: 

This  paper  deals  with  periodic  waves  on  an  ocean 
of  infinite  depth.   The  flow  is  assumed  to  be  two-dimensional, 
incompressible,  steady,  and  irrotational.   The  impossibility 
of  the  existence  of  an  asymmetric  wave  is  proved.   This  is 
accomplished  through  an  application  of  Steiner  symmetrization. 
Also  discussed  is  the  shape  of  possible  periodic  waves. 

Using  the  calculus  of  variations,  we  set  up  an 
extremal  problem  involving  the  kinetic  energy,  an  area  integral, 
and  the  potential  energy.   For  waves  of  small  amplitude  the 
kinetic  energy  is  shown  to  be  a  minimum  if  we  fix  the  area  and 
the  potential  energy.   This  is  accomplished  by  showing  the 
first  variation  to  be  zero  and  the  second  variation  to  be 
positive.   Since  the  kinetic  energy  is  closely  related  to  the 
Dirichlet  integral,  this  is  a  generalization  of  the  Dirichlet 
principle.   This  result  is  applicable  in  showing  the  existence 
of  periodic  surface  waves. 


v 


CHAPTER  1 


INTRODUCTION 


In  this  paper  we  shall  discuss  periodic  waves  at  the  free 
surface  of  an  ocean  of  infinite  depth.   The  motion  of  the  fluid 
is  assumed  to  be  two-dimensional,  incompressible,  and  irrota- 
tional.   The  only  external  force  acting  on  the  fluid  is  the 
constant  force  of  gravity. 

We  shall  use  a  moving  coordinate  system  which  travels  with 
the  same  velocity  as  the  wave.   With  respect  to  this  coordinate 
system,  the  flow  that  we  are  investigating  can  be  considered  to 
be  steady. 

Since  the  motion  is  irrotational,  a  velocity  potential 

<})(x,y)  exists.   If  u,  v  are  the  components  of  velocity  parallel 

to  the  X  and  y  coordinate  axes  respectively  at  the  point  (x,y), 

then  (f)   =  u  and  ^     =  v.   The  equation  of  continuity  (as  expressed 
X         y 

for  incompressible  fluids), 

div  (u,v)  =  0  , 

then  tells  us  that  ^   is  harmonic. 

Since  ^   is  harmonic,  we  know  (j)  has  a  harmonic  conjugate, 
which  we  call  ip.      The  curves  if/  =   constant  are  the  streamlines 
and  we  may  let  ^  =  0  represent  the  free  surface.   The  function 
C(z)  =  <|)(x,y)  +  iV'(x,y)  is  called  the  complex  potential.   The 


the  flow,  which  we  denote  by  q.  Is  then  defined  by 
the      Ion 

^-\§\  ■ 

The  shape  of  the  free  surface  is  determined  by  the  free 
boundary  condition 

where  g  is  the  acceleration  of  gravity  and  A  is  a  constant. 
This  condition  is  derived  from  Bernoulli's  law 


2 
^  +  gy  +  ^  =  const.  , 


where  p  is  the  pressure  and  p  is  the  density  at  the  point  (x,y), 
Our  free  boundary  condition  follows  by  noting  that  along  this 
curve  the  pressure  is  constant. 

To  facilitate  our  investigations  we  will  assume  that  the 
flow  has  period  one  with  respect  to  the  x-direction  and  confine 
ourselves  to  one  period  D.   The  semi-infinite  region  D  is 
assumed  to  have  vertical  side  boundaries  and  an  upper  boundary 
which  we  will  denote  by  C. 

Periodic  surface  waves  on  an  ocean  of  infinite  depth  have 
been  Investigated  by  many  mathematicians,  including  Nekrassov, 
Levl-Clvlta  and  Llchtensteln.   The  existence  of  such  waves  was 
first  proved  in  1921  by  Nekrassov  [15].   In  1925,  Levl-Civlta 
[10]  published  an  existence  and  uniqueness  proof  for  the  case 


in  which  the  quantity  ^  ^  is  sufficiently  close  to  unity, 

2Trc 
where  cd  is  the  wave  length  and  c  is  the  speed  of  the  wave. 

Levi-Civita  proved  the  convergence  of  a  power  series  in  the 
amplitude  of  the  solution.   The  first  term  of  this  series  is 
the  linearized  solution  (see  Lamb  for  a  discussion  of  the  linear 
theory  [9])«   In  1931j  Lichtenstein  used  the  theory  of  non- 
linear integral  equations  to  demonstrate  existence  and  unique- 
ness of  a  solution  if 


s-^  -  m-k  , 


27rc 

where  m  is  a  positive  integer  and  k  is  sufficiently  small. 

More  recent  results  are  described  in  the  appendix  to 
Water  Waves  by  Stoker  [l6].   Martin  has  also  contributed  to  this 
subject,  and  In  particular  in  a  paper  [3]  written  with  Dunninger 
in  1966  he  discusses  the  question  of  uniqueness. 

In  this  paper  we  shall  investigate  periodic  surface  waves 
on  an  ocean  of  infinite  depth  by  employing  a  variational 
principle.   The  first  sections  deal  with  the  construction  of 
suitable  functionals  which  we  will  need.   Specifically,  we  will 
introduce  the  kinetic  energy  associated  with  the  disturbed  part 
of  the  flow,  which  is  given  by 


=  ff   [V(V/-  cy)]^  dA  , 


M 

D 


where  c  is  the  speed  at  y  =  -00  .   We  also  consider  certain 
finite  area  and  potential  energy  integrals 


=// 


dA 
and 


// 


E  =  2g  JJ  y  dA  . 

We  will  prove  the  following  variational  formula: 

5(M  +  (?v  -c^)A  -  E)  =  /  (q^  +2gy-  A)6nds  , 

C 

which  shows  that  the  free  boundary  condition  we  obtained  from 
Bernoulli's  law  is  a  transversality  condition  for  the  extremal 

problem 

p 

M  +  (A-c)A-E  =  minimum  . 

More  precisely,  from  the  above  formula  we  see  that  if  A  and  E 
are  held  fixed,  then  M  has  a  stationary  value  when  ^  describes 
a  wave. 

The  functional  M  +(A  -c  )A  -E  is  used  to  prove  the  main 
result  of  the  paper,  which  is  the  derivation  of  a  uniqueness 
theorem  excluding  the  existence  of  certain  asymmetric  waves. 
We  will  apply  continuous  Stelner  symmetrization  (cf.  [1^])  with 
respect  to  the  y-axis  to  the  above  functional.   Prom  this  it 
will  be  shown  that  the  functional  is  convex  with  respect  to 
continuous  Steiner  symmetrization.   We  then  conclude  that  the 
wave  must  be  symmetric  with  respect  to  a  vertical  line  through 
a  trough  or  crest. 


We  will  formulate  a  more  specific  problem  in  the  calculus 
of  variations  for  f.      In  the  case  of  Levi-Civita ' s  waves  of 
sufficiently  small  finite  amplitude,  it  will  be  shown  that  the 
wave  minimizes  M  for  fixed  A  and  E.   However,  we  have  not  as 
yet  been  able  to  prove  this  result  for  the  general  case. 

Finally  we  prove  that  the  free  boundary  of  the  wave  has 
a  non-parametric  representation  y  =  y(x),  where  y(x)  is  a 
single-valued  function  of  x  with  no  vertical  tangent  at  any 
point.   To  prove  this  theorem,  we  make  use  of  Hopf ' s  lemma, 
which  gives  a  strict  inequality  for  the  normal  derivative  of  a 
subharmonlc  function  when  evaluated  at  a  boundary  point  at  which 
a  maximum  occurs. 

In  1965  Garabedlan  [5]  published  a  paper  on  periodic 
surface  waves  on  an  ocean  of  finite  depth  in  which  he  formulates 
an  extremal  problem  in  order  to  prove  the  existence  of  such  a 
wave.   Since  the  wave  described  by  ip   for  this  case  is  a  minimax 
solution  rather  than  a  minimum,  the  use  of  Morse  theory  was 
necessary. 

The  results  of  this  paper  may  be  used  as  a  basis  for  an 
existence  theorem  in  the  case  of  an  ocean  of  infinite  depth. 
In  particular,  the  result  that  the  wave  minimizes  M  for  fixed  A 
and  E  shows  that  Morse  theory  can  be  avoided.   Our  theorem  that 
the  wave  has  a  non-parametric  representation  y  =  y(x)  can  be 
used  in  restricting  the  class  of  competing  functions  involved  in 
solving  the  minimum  problem  to  those  which  have  this  same 
representation.   Symmetrization  yields  a  corresponding  result 


In  the  other  direction.   It  tells  us  that  if  we  restrict  our- 
selves to  '^he  wave,  then  the  free  surface  has  a 
non-parametric  representation  x  =  x(y)  with  Just  two  branches. 
Hence  we  may  restrict  the  class  of  competing  functions  to  those 
which  can  be  described  in  the  same  way.   Such  an  existence 
proof  by  variational  methods,  however,  is  beyond  the  scope  of 
this  paper. 


CHAPTER  2 

ESTABLISHING  FUNCTIONALS  FOR  THE  INVESTIGATION 
OF  SURFACE  WAVES 

In  this  chapter  we  look  at  the  behavior  of  the  stream 
function  ^  as  y  ->•  -00  and  then  set  up  suitable  functlonals  with 
which  to  Investigate  periodic  surface  waves  on  an  ocean  of 
infinite  depth. 

In  the  case  of  an  ocean  of  finite  depth,  Garabedlan  [5] 
used  the  Dlrlchlet  Integral  as  a  functional  for  his  investiga- 
tions of  the  existence  and  uniqueness  of  surface  waves.   Since 
the  Dlrlchlet  integral  for  the  case  of  an  Infinitely  deep  ocean 
may  diverge,  we  must  alter  it  to  obtain  a  convergent  Integral 
suitable  for  our  work.   To  do  this  we  must  first  discuss  the 
behavior  of  certain  harmonic  functions  f   defined  on  a  semi- 
infinite  region  D,  whose  upper  boundary  is  a  curve  C  and  whose 
side  boundaries  are  the  two  vertical  lines,  x  =  0  and  x  =  1. 

2.1.   Asymptotic  Behavior  at  Great  Depth 

Theorem  2.1.   Let  f   be  defined  on  a  semi-infinite  domain  D  whose 
upper  boundary  is  the  curve  C  and  which  is  bounded  on  the  sides 
by  the  vertical  lines  x  =  0  and  x  =  1  (see  Fig.  1).   We  assume 
V^  has  the  following  properties: 

t  =  0 
on  C, 

^(0,y)  :=  i^{l,Y) 
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for  all  y, 

AV'  =  0 

in  D  and  on  the  boundary  of  D,  and 

(VV^)2  =  0(1) 

as  y  -♦  -oo  . 

Through  the  techniques  of  conformal  mapping,  we  show  that 

i/     and  t(/     approach  constant  values  exponentially  fast  as  y  -»■  -oo 
X      y 

Proof:   We  may  map  the  strip  0  <  x  <  1  onto  the  w-plane 

(w  =  u  +  iv),  silt  along  the  line  v  =  0,  making  -oo  in  the  z- 

plane  correspond  to  zero  in  the  w-plane.   The  curve  C  in  the 

z-plane  (see  Fig.  1)  maps  onto  a  finite  curve  C'  in  the  w-plane 

(see  Fig.  2).   This  map  is  effected  by 


-27riz 
w  =  e 


Note  that  the  line  x  =  1  is  mapped  into  the  upper  part  of  the 
cut  and  x  =  0  is  mapped  into  the  lower  part  of  the  cut. 

Let  C(z)  =  (f)  +  1^,  where  ^   is  the  harmonic  conjugate  of  V- 
Then  C(z)  is  analytic  in  D  and  on  the  boundary  of  D.   The 
function 

may  ue  written  as 

r  (z)  =  «t)^+  i^j^ 
or 

C'(X)    =    V'y  +  lV'^    . 
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We  note  that  ^'(z)  is  also  analytic  in  D. 

Now  we  can  solve  for  z  as  an  analytic  function  of  w  in 
the  slit  w-plane.   Hence 

g(w)  =  C'(z(w))  3=  C'(z) 


I 


is  an  analytic  function  of  w  in  a  bounded  region  D  of  the  w- 
plane.  This  region  may  be  visualized  as  the  shaded  region  of 
Fig.  2. 

Let  us  look  at  what  happens  as  v  -f  0  through  positive  and 
negative  values  of  v.   Since  f   is  periodic,  we  have 

Thus 

lim  g(w)  =  lim  g(w) 


v- 


0+        V-+0- 


for  u  ^  0.      This  implies  that  g(w)  may  be  extended  to  the  points 
on  the  slit  bordering  the  domain  D'  in  the  w-plane,  except 
possibly  the  origin.   We  assumed 

lg(w)|2  =    K'(z)|  <  K^ 

near  w  =  0  or  y  =  -oo  .   Hence  g(w)  has  at  most  a  removable 
singularity  at  w  =  0.   To  make  g(w)  analytic  at  w  =  0,  we  put 


a^  =   g(0)  =  lim  g(w)  • 
"         w— 0 


Then,  it  follows  from  Taylor's  theorem  that 


CD 

g(w)  =  5~"  a^w 


in  a  neighborhood  of  w  =  0.   Therefore, 

^•(z)  =  g{w)  =  yZl   a  e-^^"^ 


5  ^"' 


for  y  sufficiently  large  and  negative.   Since  C'(z)  =  V  +  iV'  > 

y   ^ 

we  arrive  at  the  expansions 


Vy  =  Re  a^  +  Re  1^  a^ 


-27rinz 
e 


and 

GO 


V'^  =  Im  a  +  Im  rZ  a  e"^^^"^ 


Thus  we  arrive  at  the  desired  result,  i.e.. 


Vy  =  Re(a^)  +  0(e^^) 


and 


=   Im(a^)  +  OCe^"^) 


as  y  -►  -00  .   In  other  words,  the  horizontal  and  vertical  compo- 
nents of  velocity  approach  constants  exponentially  fast. 

2. II. a.   Introduction  of  the  Quantity  M 

Let  V  be  a  stream  function  defined  on  the  domain  D,  which 
is  the  same  domain  as  in  previous  sections  (see  Fig.  1).   Let 
us  assume  the  following  for  t}/: 

V'(0,y)  =  V(l,y) 


10 


for  all  y, 


IV^I^  =   0(1) 


as  y  -►  -00  ,  and 

^  =  0  on  the  free  surface  . 

Then  it  follows  from  Theorem  2.1,  that  f   has  the  following 
expansion  near  y  =  -oo  : 

t  =  cy  +  y^  +  >    a^^e  '^^^   cos  (27rnx+  \>^)    . 
n=l 

Note  that  the  a  here  are  not  quite  the  same  as  the  a  used 
n  ^  n 

previously. 

Since  we  have  investigated  the  behavior  of  f   near  y  =  -od  , 
we  can  now  introduce  the  quantity  M.   This  quantity  M  will 
represent  the  kinetic  energy  associated  with  the  disturbed  part 
of  the  flow  and  will  be  an  integral  involving  the  term  (V^) 
plus  other  terms  which  are  added  on  to  make  it  convergent  over 
the  domain  D.   More  precisely,  we  define 


=  JJ    [S/{f  -cy)f   dA  . 


M 

D 


It  is  easily  seen  that  this  integral  is  convergent,  since 
we  have  investigated  the  behavior  of  tl/   and  (V^)   near  y  =  -oo  . 
After  subtracting  the  term  cy  from  Tp ,    we  see  that  the  square  of 
the  gradient  of  the  resulting  function  dies  out  exponentially 
as  y  -♦  -00  . 


11 


Most  of  the  results  In  this  paper  rest  upon  the  investl- 
gatl      properties     .   Formulas  for  the  first  and  second 
varlatiw.io  ^f  M  will  be  derived.   The  functional  M  and  its 
first  variation  6M  will  be  important  in  showing  the  Impossibility 
of  the  existence  of  certain  asymmetric  periodic  waves.   In  the 
special  case  in  which  the  wave  is  almost  flat,  we  indicate  in 
Chapter  4  how  M  might  be  used  to  set  up  an  existence  proof  for 
periodic  waves  on  an  ocean  of  infinite  depth.   In  particular  we 
show  that  for  a  fixed  area  and  a  certain  fixed  potential  energy, 
the  second  variation  of  M  is  positive. 

2.II.b.  The  First  Variation  of  M 


Suppose 

J(e)  =^F(x,y,^,^^,^y)dA  , 


D 


where  the  domain  D  and  the  functions  ^,  rj/   ,   and  ^  depend  on 
the  parameter  e,  and  e  is  small.   Then  the  first  variation  of 
the  functional  J(e)  can  be  defined  by  the  formula 


de 


e 
e=0 


Let  D  be  the  same  semi-infinite  domain  as  mentioned  in 
previous  sections.   Its  upper  boundary  is  C;  its  side  boundaries 
are  the  vertical  lines  x  =  0  and  x  =  1.   We  obtain  the  domain 
D*  from  D  by  varying  only  the  upper  boundary  of  D.   The  upper 


12 


boundary  C*  of  D   is  obtained  by  making  a  variation  of  the  form 

5n  =  ep(s) 

along  the  inner  normal  to  C  with  respect  to  the  domain  D.  The 
function  p(s)  is  assumed  to  be  a  real  analytic  function  of  the 
arc  length  parameter  s  on  C. 

It  is  assumed  that  ^  is  a  symmetric  stream  function,  that 
is. 


in  D, 
on  C, 

for  all  y,  and 

as  y  — »•  -00  . 


A^  =  0 


f  =   0 


^(0,y)  =  ^(l,y) 


(V^)2  =  0(1) 


Let  Tl/     be  the  harmonic  function  defined  on  D*  which 
satisfies  the  analogous  conditions  with  respect  to  D*,  that  is, 


* 


In  D*, 
on  C*, 

for  all  y,  and 


A^/  =  0 


f*  =   0 


1//  (0,y)  =  ^*(l,y) 


15 


(V^*)2  =  0(1) 
as  y  -♦  -OD  . 

.._  :.--■_  that  the  curve  C  is  assumed  to  be  analytic  and 
that  the  variation  is  analytic.   This  implies  that  C  is 
analytic.   Hence  if  C  and  C*  are  sufficiently  close,  then  by 
Schwarz's  principle  of  reflection  V  Sind  ip     can  both  be  extended 
outside  of  their  domains,  if  necessary,  so  that  they  are  both 
harmonic  in  D  and  D  and  on  C  and  C*. 

Suppose,  at  first,  we  make  the  assumption  that  D*  lies 
entirely  in  the  interior  of  D.   Let  n  be  the  inner  normal  to  D, 
and  let  D,  be 

D  n  (y  i  y  >  -k,  k  >  0)  . 
Then  we  may  rewrite  M  as 


M  =  lim  //  [V(^-  cy)]^  dA  . 


We  have  through  an  application  of  Green's  theorem 

=  J   cy  "5^  (^  -  Gy)d3  , 


M 


since  the  line  integrals  over  the  vertical  boundaries  of  D, 
cancel  each  other  and  since,  in  the  limit  as  k  -♦  oo  ,  the  line 
integral  over  y  =  -k  vanishes. 

Let  A  be  the  area  bounded  by  the  curve  C  and  the  x-axls, 
Then  M  can  be  written  as 


1^ 


M  =  c  r  y  II  ds  +  c^A  , 


and  therefore 


M  -  c 


(^l^-*i'^^ 


^  C 


2 


This  follows  because  ip  =  0  on   the  curve  C  and  hence  the  second 
Integral  on  the  right  of  the  above  equation  is  zero.   Now 

//   V^-V^*dA  =  /  (V^)^  5nds  +  O(e^)  . 
D-D*  C 

On  the  other  hand,  by  Green's  theorem. 


Vi//-VV/*dA 


^     an 


11/     -^  ds 
^      on 


D-D 


The  second  integral  on  the  right  side  of  the  equation  vanishes 

on   C  .   Adding   /  f  -^^   ds  to  the  first  integral  on  the  right 

C 
hand  side  of  the  equation,  which  leaves  it  unchanged  because 

^  =  0  on  C,  we  get 

ff    VV/.V^-dA  =  -/('/'*  II  -  ^||^)ds 


D-D^ 


^      an        ^   an 


y=-k 


Now 


00 


r  =  yI  +  cy  +Y~  ''n^  c°s  (n7nc  +  ^*) 


and 


n=T 

00 


V'  =  y^  +  cy  +  S        a  e  ^^  cos  (n-rrx  +  h    ) 
n=l 
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near  y  =  -oo 


substituting  for  rj/,   ^,  rj/* ,   and  ■^,    and 
.3  of  the  expansion,  we  have 


D-D* 


VV'«V^*dA 


J. 


C-[y*+cy]c   +   [yQ  +  cy]c)dx 


0 
y=-k 


=  --(y^yo^  • 


We  now  derive  the  formula 


M  -  c  A  =  -cy^  , 


which  relates  the  kinetic  energy  M  with  the  coefficient  y^  in 
the  Taylor  expansion  for  the  stream  function  tJ/   near  y  =  -oo  .   An 
application  of  this  formula  allows  us  to  get  the  variational 
formula. 

We  have 


M  =  c 


f 


^  C  C 


■Tp   1^  ds 
^  dn 


2 
+  c  A  . 


Thi--.  ■—'"--  Green's  theorem  and  taking  into  account  that  y  and 
Tp   are  harmonic  in  D,  we  obtain 


M  =  c 


(-y|i.*|i)dx 


0 
y=-k 


2 

+  c  A 


Substltutlnc   the  above   expansion   for  ip  near  y  =   -oo  ,    I.e.,    for 
k   suri'icientiy   large,    give:; 


16 


M  =  c 


± 
f    (  (-y)(-c)  -  (cy +  y^))dx 


0 
y— k 


2 

+  c  A  , 


which  reduces  to 


or 


M  =  -cy  +c  A  , 


M  -  c  A  =  -cy   . 


Since  we  may  write  the  analogous  formula 


M*-  c^A*  = 


-cy 


o  ' 


the  total  variation  of  M  -  c  A  is 


Also,  we  have 


(M*-M)  -c^(A*-A)  =  -c(y^-  y^) 


M* -M  =  -c(y*-  y)  +c^(A*-  A)  . 


From  the  equations 


and 


we  have 


Therefore 


//   V^-V^*dA  =  /  (V^)^  5nds  +  O(e^) 


D-D' 


JJ    V^.V^^dA  =  -c(y*  -y^)  , 


D-D 


-c(y*-  y^)  =  j  (V^)^  5nds  +  0(£^) 
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(M*  -M)  -  c^(A*  -A     ^  e^) 


C 

rq^(5n)ds  +  0(e^; 


C 


where  q  is  the  speed  along  C.   Denoting  the  first  variation  by 
the  symbol  5,  and  noting  that 


6A  =  -  /  (5n)ds  , 


we  obtain  the  final  result 

5M  =  '  (q  -  c  )5nds 


We  have  assumed  that  the  domain  D  is  interior  to  D. 
However,  this  is  not  necessary.   If  one  of  the  domains  is  not 
inside  the  other,  we  consider  a  third  domain  D   which  has  the 
same  vertical  side  boundaries  as  D  and  D  but  is  contained 

vi  oth  D  and  D  .   Let  C   denote  the  upper  boundary  of  D 

**   * 
We  have,  by  the  above  computation,  expressions  for  M   -  M  and 

M   -  M  up  to  first  order.   If  we  denote  the  variation  along  the 


inner  normal  to  C  by  the  symbol  5n  and  the  variation  along  the 
inner  normal  to  C  by  the  symbol  5n,  where  6n  and  5n  both  bring 
us  to  a  point  on  the  curve  C   ,  then  we  can  write  5n  as 


5n  =  6n  -  5n  +  O(e^) 


Up  to  first  order  in  e,  we  have 
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•v 


and 


M**-  M  =  /  (q^-  c^)5nds 


We  can  rewrite  the  second  equation  in  the  form 

M**  _M*  =  /  (q^  -  c^)5nds  +  O(e^)  . 

C 

Then 

M*  -  M  =  (M**-M)  -  (M**-M*) 

=  j  (q^  -  c^)5nds  +  O(e^)  . 
C 

Hence  the  first  variation  of  M  is  given  by 
Theorem  2. II. 

5M  =  /  (q^  -  c^)6nds  . 


2.II.C.  The  First  Variation  of  M  for  a  Wave 

In  Section  2.II.b.,  we  derived  the  formula 

m  =         (q^  -  c^)5nds 
C 

for  the  first  variation.   To  derive  this  formula,  we  made  use 
of  the  functional 


A  =JJx(x,y)dA  , 


where 
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X(x,y)  =  1 

for  y  >■  -  awU 

X(x,y)  =  0 

for  y  ^  0,  and  the  first  variation  of  A 

5A  =  -  /  6nds  . 
C 

We  now  assume  that  Tp   describes  a  wave  and  define  the 
functional 


where 


P(e)  =  M  +  (X  -  c^)A  -  E  , 


E  =  JJ  2gY   x(x,y)dA  . 


The  integrals  M,  A,  and  E  are  defined  over  the  domain  D  ,  which 
depends  upon  e.   This  domain  is  described  in  detail  in  Section 
2.II.b.   The  i'-'^^'^-^al  E  represents  the  potential  energy  of  the 
flow.   The  first  variation  of  E  is  obviously  given  by 

6E  =  -  /  2gy  6nds  . 
C 

Since  C  is  the  free  boundary  of  the  wave,  we  no  longer  have  to 
require  that  C  be  analytic  in  order  to  derive  the  variational 
formulas,  because  the  requirement  can  be  deduced  from  the  free 
boundary  condition  [11] 

q^  +  2gy  =  A  . 
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Since  we  know  &A,  5M,  and  5E,  we  may  now  write  down  an 
expression  for  the  first  variation  of  F.   This  Is  given  by  the 
formula 


=  J  (q^+  2gy  -  A)6nds 


6F 

C 

Therefore  a  necessary  condition  for  F(£)  to  have  a  stationary 
point  when  e  =  0  is  that 


q^  +  2gy  =  A  , 


which  is  the  free  boundary  condition.   In  other  words  ip   satisfies 
the  free  boundary  condition  and  hence,  describes  a  wave.   This 
result  is  the  basis  for  the  uniqueness  theorem  of  the  next 
chapter  and  for  the  setting  up  of  the  variational  principle  of 
the  concluding  chapter. 
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CHAPTER  3 
THE  SHAPE  OF  PERIODIC  WAVES 

This  chapter  deals  with  the  possible  shapes  of  periodic 
waves.   If  we  assume  that  all  the  level  curves  V(x,y)  =  constant 
rise  and  fall  just  once  in  one  period,  it  will  follow,  from  an 
application  of  Steiner  syrametrization,  that  no  asymmetric 
periodic  wave  is  possible.   This  is  the  main  theorem  of  this 
section. 

By  assuming  that  the  free  surface  of  the  wave  ^  =  0  rises 
and  falls  just  once  in  one  period  we  can  show  that  the  curves 
V'(x,y)  =  t  have  the  same  property  providing  t  is  sufficiently 
large  and  negative.   Thus  the  hypothesis  of  the  theorem  may  be 
weakened.   With  an  additional  restriction  on  the  free  surface 
alone  we  may  also  prove  that  all  curves  ^  =  t  rise  and  fall 
precisely  once  in  one  period. 

It  will  also  be  shown  that  a  periodic  surface  wave  can  be 
described  non-parametrically  by  a  single-valued  function  y  =  y(x) 
which  has  no  vertical  tangent.   This  last  result  will  depei 
Hopf's  lemma  for  subharmonic  functions,  which  deals  with  the 
normal  derivative  at  a  point  on  the  boundary  of  a  domain  at  which 
a  maximum  occurs  [1].   These  results  have  an  application      e 
study  of  existence  and  uniqueness  of  surface  waves  by  variational 
methods. 
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3.1.  The  Non-Exlstence  of  Certain  Periodic  Asymmetric  Waves 

We  proceed  to  prove  the  non-existence  of  periodic 
asymmetric  waves  under  the  assumption  that  all  the  streamlines 
ascend  and  descend  just  once  In  a  single  period.   In  other 
words,  we  assume  that  each  streamline  ^(x,y)  =  constant  has  a 
non-parametric  representation  x  =  x.(y)  in  the  domain  D 
describing  one  period  of  the  flow,  with  just  two  branches  x, 
and  Xp. 

Motivated  by  Section  2.II.C.,  we  shall  set  up  a  functional 
F  which  will  depend  on  a  parameter  9,    where  0  ^  0   <   2.      In  this 
interval  F{6)   will  be  shown  to  be  convex  and  the  derivative  of 
F(0)  will  be  shown  to  vanish  for  9=0  and  9=2.      Thus  it  will 
be  proven  that  F(9)  must  be  constant.   This  conclusion  will 
result  in  a  contradiction  unless  the  streamlines  ^  =  constant 
are  symmetric  with  respect  to  a  vertical  line  through  a  crest. 

We  recall  that  D  is  the  semi-infinite  domain  which  lies 
between  x  =  0  and  x  =  1  as  in  previous  sections.   In  this  section, 
the  function  ^  defined  on  D  is  assumed  to  describe  a  wave.   For 
now,  we  translate  the  y-axis  to  the  right  by  l/2  (see  Fig.  3). 

Suppose  ■4/{x,y)    is  asymmetric  with  respect  to  the  new 
y-axis.   Then  ?//(x,y)  and  ^(-x,y)  define  different  streamlines. 
This  means  that 

x-[(y)  f   -X2(y)  . 

We  proceed  by  constructing  a  new  set  of  streamlines  with  two 
branches  x  =  X^(y)  and  x  =  X2(y)  defined  by 


(1) 


H-^.-<'-K^ 


/x.  +x„ 


Where  0  <   9  <  2. 

For  each  value  of  0,  we  may  define  a  function  V(^)  whose 
level  curves  have  the  branches  X,  and  X^  corresponding  to  the 
same  value  of  9.      The  function  Tp{6)    is  defined  on  the  domain  Dg 
with  the  same  vertical  side  boundaries  as  D  and  with  the  curve 
^(0)  =  0  for  -1/2  j^  X  <  1/2  as  upper  boundary.   Note  that  i/{d) 
is  not  necessarily  harmonic.   We  may  Interpret  the  transforma- 
tion (1)  as  continuous  Steiner  symmetrization.   For  9=0,  we 
have  the  identity  transformation.   For  9=1,  the  transforma- 
tion is  Steiner  symmetrization  with  respect  to  the  y-axis.   For 
0=2,  ^{x,y)  is  mapped  into  V'(-x,y). 

Let 

F(e)  =  M  +  (A-c^)A  -E 

=  ff[[^{f{d)  -cy)f  +   {A-c2)x(x,y)  -  2gy  x(x,y)]dA  , 
^9 


where,  again. 


for  y  >•  0,  and 


X(x,y)  =  1 
X(x,y)  =  0 


:  :•  y  <  0.   This  functional  F(0)  is  well-defined  provided  we 

^ hat  M(0)  converges  for  each  9,    where  0  <   9  <  2.      The 


2k 


convergence  of  M(9)  will  be  shown  trivially  once  we  have  shown 
it  to  be  convex. 

The  crucial  part  of  this  proof  consists  in  showing  that 
F(9)  is  convex.   It  will  be  convenient  to  use  the  following 
formulation  of  convexity:   F(0)  is  convex  if 


Y!l^^  <  I  [f(0-l)  +F{e^)]    . 


The  functionals  A  and  E  are  obviously  independent  of  0,    because 
the  transformation  (1)  is  simply  a  horizontal  shift  of  the 
curves.   Therefore,  to  prove  the  convexity  of  F,  it  is  sufficient 
to  show  M(0)  is  convex.   To  accomplish  this,  we  rewrite  M  as  an 
Iterated  integral  and  show  that  the  inner  integral  with  respect 
to  X  is  a  convex  function  of  Q. 
We  note  that 


M(0)  =  f  [f   [^(^-  cy)]^dx 


dy 


Let  us  now  define  the  functional 

m(0,y)  ^  f  [V(^(0)  -  cy)]^  dx  , 

where  m(0,y)  forms  the  integrand  of  M(0),  i.e., 

i(e)  =  J  m{e,Y)dy    . 


Ml 


We  proceed  now  to  establish  the  convexity  of  m(e,y).   The 
level  curves  of  ip   are  the  curves  i//(x,y)  =  t,  where  the  quantity 
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t  varies  from  zer.     -  >  .   We  make  a  change  of  variables  by 
introducing  t  as  :.     -.ependent  variable  replacing  x.   The 
variable  x  will  then  depend  on  t  and  y,  thus  becoming  a 
dependent  variable. 
Then  we  have 


V'^(0)^dx  =  -^^   dt 


^ 


if  we  are  on  the  Jth  branch  of  the  curve  V(^)  =  t.   Note  J  =  1,2. 
Similarly, 


2^     2 
c  dx  =  c 


bx 


'St 


J 


dt 


and 


^je)^dx=^(^rdt  . 


^F 


Expanding  the  expression  for  m(9,y)   gives 


m(0,y)   =  f  [i^^{e)^  +  ^^{d)^  -2cV/y(0)  +c^]dx 


Thus,  by  making  the  above  substitutions  and  letting  a  prime 
denote  partial  differentiation  with  respect  to  y  and  a  dot 
denote  partial  differentiation  with  respect  to  t,  we  get 


m(e,y)  =  YZ. 


^(x;)^+i     2,-  , 

— J-^ +  c^|X  I  -  cX;  sgn  X 


dt 


We  now  determine  the  signs  of  X^  and  X^.   Note  that 
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and 


*i  -  K^¥^^ 


^2  -  ^2^^^ 


The  quantity  4|-  is  always  positive  since  t  varies  from  0  to  -co 
as  we  descend  from  the  upper  boundary  of  D  to  y  =  -oo  .   Now, 
recall  that 

and  (1) 


where  0   <   Q   <   2.      The  branch  x-,  (y)  of  ip   =  t  was  chosen  so  that 
it  always  has  positive  slope  and  the  branch  x^Cy)  was  chosen 
so  that  it  always  has  negative  slope;  that  is 


xj(y)  >  0   and   x^(y)  <   0  . 


We  see,  by  differentiation  of  equations  (1)  with  respect  to  y, 
that 


and 


X-i  —  x-j  -  -p  vXq  +  yir^ ) 


■X-p  —  Xp  -  p-  ^X,+Xpj 


where  0   <   9   <^  2.      Hence  we  may  write 

X-,  =  (1-  ■^■^^x  ~  '^   ^2 

and 

I  _  ,    e>  I    0   1 
Xp  —  \  i  -  '^'■^o    ~    2  ■'^l  ' 
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a 

Where  0  <  1  -  -^  Jl  1 
This  implies  that 


and 


Thus  it  is  obvious  that  X '  >  0  and  Xp  <  0. 


X^  >  0 


X^  ^  0 


We  rewrite  m  by  substituting  +1  for  sgn  X,  and  -1  for 
sgn  Xp.   We  also  rewrite  |X^|  as  X,  and  |Xp|  as  -Xp.   Thus  we 
have 


m(0,y)  =  f 


(X,  )^+l    p.      , 

— —. +  c'^X,  -  cX, 

X.        ^    -^ 


dt 


/ 


^(X')^  +  1    p. 
—^ c^'Xg  +  cXg 


-X, 


dt  . 


We  now  proceed  to  show  the  convexity  of  m  with  respect  to 
6   by  making  use  of  the  linearity  of  the  transformation  (1)  in  9. 
Let 


^1,1  =  ^l(^)  '        ^1,2  =^1(^2) 


Then 


and 


^2,1  =  ^2^^1^  '        ^2,2  =  X2^^2^  • 


^1  I — 5 )=  ^   ^^1,1  "^  ^1,2^ 


Xgl  -^-^y  ^  ^^2,1  +  ^2,2^ 


since  we  have  linearity  in  9.   Substitution  in  the  equation  for 
m(0)  gives 


m 


h^'2 


yj  =LZ  (-1) 


j+i 


X".  ,  +  X'.  ^\  2 


+  1 


J>1     J, 2 


+  c 


2^^j,l  +  ^j,2 


-  c 


^Li'^^.U 


dt 


We  may  rewrite  the  non-linear  terms  In  the  above  Integral  as 


(-1) 


J+1 


X'.  ,  +  X'.  ^N2 


+  1 


J>1     J, 2 


-  (-1) 


,H-iCii.i!!ii,i^iZ^ 


2X.  , 


2X 


J, 2 


♦      •     p     .     I       1   .     p 

(X.  T  -X.  ^)   +  (X.  ,x.  o  -  X.  tX.  o) 
J,l    3,2  J,l  J, 2     j,l  J, 2^ 


2X.  ,X.  ^(X.  .  +  X.  o) 
j,l  J, 2'  j,l     J, 2' 


This  is  done  by  completing  the  square.   Since  the  other  terms 
in  the  integrand  are  linear  in  9,  we  have 


(2)   m  f-l^,y]=  \   m(0^,y)  +  \  m{d^,Y) 


-LZ(-i) 


J+1 


(X.  -,  -X.  o)^+  (X.  .x'.  o-x'.  .X.  ^y 

J>1    3,^  J,l  J, 2    j,l  J, 2^ 


2^J,lXj,2('^J,l+^J,2> 


dt 


Thus 


(5) 


m(  ^  g  ^  )  <  I  [m(e^,y)  +  m(02,y)]  > 
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which  means  that  m(e,y)  Is  convex.   The  method  used  in  the 
^^...-.-^..  --  -..is  is  equivalent  to  an  application  of  Schwarz's 
inequality. 

We  are  now  ready  to  prove  that  F(0)  is  well-defined  by 
showing  the  convergence  of  M(0).   We  note  that 

m(0,y)  >   0  , 

since  its  integrand  is  always  positive.   From  the  convexity  of 
ni(0,y),  we  may  conclude  that 

in(0,y)  <  max  (m(0,y )  ,m(2,y ) ) 

<  m(0,y)  +  m(2,y)  . 


We  recall  that 


M(0)  =Jm{9,y)dy    , 


and  therefore 


^0    +e    N       n     /01+02      \ 

which  means  that  M(0)  is  convex.   Hence,  we  have  established 
the  convexity  of  F(0).   Since  M(0)  =  /  m(0,y)dy,  we  have 

0  <  M(0)  <  M(0)  +  M(2)  . 
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Because  M(0)  and  M(2)  converge,  we  conclude  that  M{6)   must 
converge  for  0  ^  9   <_  2.      Thus  the  functional 

F(e)  =  M  +  (A-c^)A  -  E 

Is  well-defined. 

We  may  separate  the  formula   for  6F,  the  first  variation 
of  F,  into  two  parts,  a  part  arising  from  the  variation  5?//  of 
the  stream  function  ip   over  the  fixed  domain  D  and  a  part  corre- 
sponding to  the  normal  variation  6n  of  the  boundary  of  the 
domain  D.   Thus 

5F  =  /y  5[V(^  -  cy)]^  dA  +  /  (q^  +  2gy  -  A)5nds  , 
D  C 

where  the  second  term  of  the  formula,  representing  the  part  from 
varying  the  domain,  was  derived  in  Chapter  2.   For  9=0   and 
0=2,  the  first  term  vanishes  through  an  application  of 
Dlrichlet's  principle,  while  the  second  term  vanishes  because 
^(x,y)  and  ^(-x,y)  satisfy  Bernoulli's  law 


q^+  2gy  -  A 


on  the  free  boundary,  with  the  same  constant  A.   Thus 

F'(0)  =  F'(2)  =  0  . 

We  may  conclude,  therefore,  that  F(0)  Is  a  convex  function 
whose  derivative  F'(9)  vanishes  at  the  end  points  of  the  interval 
0   <   9   <   2.      Hence  F(0)  must  be  constant  in  this  interval. 
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Returning  t  ,  we  see  that  the  F(0)  bel:. 

for  j  =  1,2  ex. .a  I'or  all  d,    where  0  <   9  <  2.      I:.  ^-.  .  .  ..^r 


However, 


and 


Thus 


X^(0)  =  X^{2)  . 


X^(0)  =  x^ 


X^(2)  =  -Xg  . 


x^(y)  =  -X2(y) 


This  gives  us  a  contradiction  unless  ^(x,y)  =  t  and  ^(-x,y)  =  t 
define  the  same  curves  for  all  t.   Hence  the  flow  V  has  symmetry 
with  respect  to  a  vertical  line  through  a  crest,  or,  equiva- 
lent ly,  through  a  trough. 

3. II. a.  A  Theorem  on  the  Behavior  of  Streamlines  at  Great  Depth 

Theorem  3. II. a.   Let  V'(x,y)  be  a  periodic  stream  function  having 
period  one  with  respect  to  its  x-coordinate  and  such  that 
[VTp)      =   0(1)  as  y -♦  -QD  .   Assume  the  curve  Tp{x,y)   =  0   ascends 
and  descends  Just  once  in  a  period.   Then  the  level  curves 
V'(x,y)  =  t  have  the  same  property  for  t  large  and  negative. 

Proof;   Let  ^^   represent  one  period  of  the  wave  (see  Fig.  4). 
We  note  that  the  side  boundaries  of  Jb   are  not  necessarily 
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vertical  as  for  the  domain  D  described  In  previous  sections. 
The  upper  boundary  of  J^  ,  which  we  denote  by  C,  will  fall  on 
the  level  curve  ^  =  0. 

Using  the  complex  potential  t,,    we  map  the  domain  J^   con- 
formally  onto  the  seml-lnf  Inlte  rectangle  t//<0,  0<(t)<cin 
the  C-plane  (see  Fig.  4),  where  c  (the  speed  at  y  =  -co  ) , 
represents  the  amount  <t)(x,y)  Increases  as  we  traverse  one 
period  of  the  wave. 

We  map  the  rectangle  plus  the  part  of  the  boundary  falling 
on  the  curves  t^  =  0  and  (j)  =  0  onto  the  unit  disc  in  the  C-^- 
plane  by  the  transformation 

-iC—    (^-i<t))— 
C^  =  e      =  e 

(see  Fig.  5).   The  point  t//  =  -qd  corresponds  to  the  origin  in 
the  C-,  -plane  and  the  curve  Tp  =  0   for  0  _^  (()  <  c  maps  one-to-one 
onto  the  unit  circle  |C-,|  =1.   Let 


0.  =  qe-^^  .^  (^  +  1^)  . 


The  angle  0   can  then  be  interpreted  as  the  angle  of  Inclination 
of  the  velocity  vector  with  respect  to  the  positive  x-axls. 

We  consider  the  image  of  the  unit  disc  U^|  <  1,  in  the 
plane  of  log  go  (see  Fig.  5).   Without  loss  of  generality,  we 
make  the  assumption  that  the  speed  at  y  =  -oo  is  one.   Then  the 
origin  in  the  C^-Plane  maps  onto  the  origin  in  the  plane  of 
log  0).   Let  us  call  the  image  of  the  unit  circle  U-j_|  =  1,  C". 
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We  use  the  free  surface  condition. 


q^  +  2gy  =  X  , 


and  the  assumption  that  C  rises  and  falls  exactly  once  tc  ..,^v; 
that  the  mapping  of  C  onto  C"  is  one-to-one.   For  9  >■  0,  y  is 
monotonically  increasing  as  we  traverse  C,  which  implies  q  is 
monotonically  decreasing.   Similarly,  for  0  <  0,  y  is  mono- 
tonically decreasing  as  we  continue  along  C,  which  implies  q 
is  monotonically  increasing.   Thus  we  may  write 

9   =   e(log  q)  , 

where  this  function  has  exactly  two  branches;  one  for  d  >  0   and 
one  for  9   <  o.      This  means  that  the  mapping  of  C  onto  C"  is 
one-to-one.   Since  the  inverse  mapping  of  lC-,|  =  1  onto  C  is 
one-to-one,  the  map  of  1C]_|  =1  onto  C"  is  one-to-one. 

Now  we  may  apply  the  argument  principle  to  conclude  that 
any  point  log  oo  ,  inside  the  domain  bounded  by  C"  has  exactly 
one  pre-image  in  the  interior  of  the  unit  circle  in  the  ^^- 
plane.   The  pre-image  of  the  real  axis  segments  OA,  OB  in  the 
plane  of  log  cu  are  two  curves  OA,  OB  In  the  /^^^-plane  inter- 
secting at  the  origin. 

We  now  show,  for  t  sufficiently  large  and  negative,  that 
the  image  of  a  level  curve  ^  =  t  in  the  plane  of  log  to  inter- 
sects 0  =  0  exactly  twice.   Since  cu  is  analytic  near  C-,  =  0,  we 
may  write 

n=l 


5^ 


In  a  neighborhood  of  C-,  =  0.   We  note  that  a.,  cannot  be  zero 
because  in  that  case  the  mapping  would  not  be  one-to-one  in  a 
neighborhood  of  C-i  =  0.   Since  conformal  maps  take  small  circles 
into  almost  circular  curves,  for  t  large  and  negative,  we  may 
say  that  the  image  of  a  circle  ip  ^   t   In   the  C-,-plane  will  be  an 
almost  circular  curve  in  the  plane  of  log  (u.   Hence  the  image 
oT  f  =  t   In   the  plane  of  log  co  will  cut  the  line  0=0  exactly 
twice.   This  means  the  image  curve  can  be  divided  into  two  arcs, 
one  on  which  6   is  positive  and  one  on  which  9   is  negative.   We 
may  now  conclude  that  a  level  curve  f  =   t   rises  and  falls 
exactly  once  in  a  period  for  t  large  and  negative.   This  is  what 
we  set  out  to  prove. 

3.II.b.  A  Discussion  on  the  Behavior  of  Streamlines  of  ^ 

In  this  section  we  make  an  additional  restriction  on  the 
upper  boundary  C  of  the  domain  G,  which  represents  one  period 
of  the  flow.   This  will  enable  us  to  obtain  the  result  of  the 
previous  section  for  all  the  level  curves  of  ip ,    that  is,  if  C 
rises  and  falls  once,  then  for  all  t,  V  =  t  rises  and  falls  just 
once  in  one  period.   The  additional  restriction  that  we  place 
on  C  is  that  its  image  in  the  plane  of  log  o),  which  we  call  C", 
bounds  a  starlike  region  with  respect  to  the  origin  of  that 

plane. 

Since  the  region  that  C"  bounds  is  starlike,  the  inner 
product  of  the  normal  to  the  curve  C"  and  the  radius  vector  is 
always  >   0.   Since  ijj  =   0  on   C"  and  (log  q,9)    Is  a  radius  vector, 
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we  see  that 


<i°« '' T^lr^  - «  i  ^  ° 


Log  q 

We  note  that  ^^  is  a  Green's  function  on  the  domain  bounded  by 
C"  and  that 

is  harmonic  in  this  domain.   The  latter  holds  because  h  is  the 
imaginary  part  of  an  analytic  function  In  this  domain,  except 
possibly  at  the  origin.   We  note,  however,  that  this  singularity 
is  removable. 

Since  h  ^  0  on  C",  it  follows  from  the  maximum  principle 
that  h  >  0  inside.   Thus,  for  any  t,  the  image  of  a  level  curve 
V'  =  t  is  starlike  with  respect  to  the  origin  in  the  plane  of 
log  CD.   This  means  that  any  line  through  the  origin  in  the  plane 
of  log  0)  will  intersect  the  image  of  ^  =  t  twice.   In  particular, 
9        "  cuts  the  Image  twice. 

We  conclude  that  ^  =  t  can  be  divldied  into  two  segments, 
cne  on  which  0   is  positive  and  one  on  which  6   is  negative.   This 
means  %/  -   z   rises  and  falls  once  in  one  period,  which  is  the 
desired  result. 

""    .   A  Theorem  on  the  Non-Parametric  Representation  of  the 
Free  Streamline 
Thrcu^rh  an  application  of  Hopf's  lemma  for  subharmonlc 
..oii£  ii],  we  will  derive  the  theorem  that  the  free  stream- 
line •-     has  a  non-parametric  representation  y  =  y(x),  where 
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y(x)  is  a  single-valued  function  of  x  with  no  vertical  tangent 
at  any  point.   This  means  the  wave  cannot  bend  back  on  itself 
and  thus  is  convex  with  respect  to  the  y-direction. 

In  the  next  chapter  we  will  discuss  a  variational  principle 
which  asserts  that,  for  fixed  A  and  E,  M  is  a  minimum  when  ij/ 
describes  a  wave.   If  we  restrict  ourselves  to  one  period  of  the 
wave,  then  by  symmetrization  we  can  prove  that  the  free  surface 
of  the  solution  to  the  minimum  problem  should  have  a  non- 
parametric  representation  x  =  x.(y)  with  just  two  branches  which 
are  symmetric  with  respect  to  the  y-axis.   These  two  theorems 
on  the  non-parametric  representation  of  the  free  surface  of  the 
wave  enable  us  to  restrict  the  class  of  competing  functions  for 
the  variational  principle  and  thus  contribute  to  an  existence 
proof  for  the  wave  problem  by  variational  methods. 

We  proceed  to  the  proof  of  the  theorem.   As  in  previous 
sections,  the  domain  D  represents  one  period  of  the  flow 
described  by  the  stream  function  ^(x,y).   For  this  theorem,  we 

assume  ^(x,y)  describes  a  wave  and  hence  satisfies  Bernoulli's 

2 
law,  q  + 2gy  =  A,  on  ^  =  0,  the  upper  boundary  of  D. 


Let 


f(x,y)  =  q^  +  2gy-  A 


2     2    2 

throughout  the  domain  D,  where  q  =  f    +  if/   .      We  have 

^   y 


f (x,y)  -*  -00 


2 

as  y  -►  -GO  ,  since  q  is  bounded  and  A  is  a  constant.   Also, 

f  (x,y)  s.  0  along  the  free  boundary  C,  and 
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Af  >  0 

In  the  domain  D;  in  other  words,  f  Is  subharmonlc.   Hence,  by 
the  maximum  principle  for  subharmonic  functions, 

f  <  0 

throughout  D. 

Let  n  denote  the  inner  normal  to  the  free  boundary  C  and 
let  9   denote  the  angle  that  the  tangent  to  C  makes  with  the 
positive  X-axis.   We  let  s  denote  the  arc  length  parameter 
along  C,  where  s  =  0  at  the  left  hand  end  point  of  C,  and  define 
the  curvature  n   by  the  formula 

d0 
ds 

Differentiation  of  f  with  respect  to  n  gives  us 

|f  =  2q  |a  +  2g  |y   . 
on    ^  on  ^   dn 


Using  the  relation 


3l  =  -"^  . 


we  have 


i  =  -^^''  ^  ^8  Is  . 


We  now  state  Hopf ' s  lemma,  which  we  will  apply  to  the  function 

Sf 

3—:   Let  D  be  a  domain  bounded  by  a  smooth  curve  and  let  u  be 

on 

a  non-positive  subharmonic  function  defined  on  D,  and  having 
continuous  first  partial  derivatives  on  the  boundary  of  D.  Let 
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n  denote  the  Inner  normal  along  the  boundary  of  D.   If  u  =  0 
at  a  point  on  the  boundary  of  D,  then  y^  Jf.  0  at  that  point. 
Unless  u  =  0,  we  get  strict  inequality.   Applying  Hopf ' s  lemma, 
we  obtain 

|I<  0  . 
an 


Hence 


2g  IJ  <  ^q^ 


This  means  that  at  any  point  along  C  where  the  tangent  is  verti- 
cal, the  curvature  must  be  strictly  positive. 

Suppose  the  curve  C  has  a  vertical  tangent  at  some  point 
P^  along  it.   Then  without  loss  of  generality,  we  may  take  the 

normal  at  P,  in  the  positive  x-direction.   Hence  -^  =  0.   This 
-L  on 

Implies,  via  the  above  inequality,  that  K   must  be  positive 
at  P^.   Since  k    Is  continuous  along  C,  fc   must  be  positive  in 
a  neighborhood  of  P^.   The  only  way  for  this  to  happen  is  for 
the  curve  C  to  bend  back  on  itself,  i.e.,  it  must  be  that  if 
x^  is  the  x-coordinate  of  P^,  then  there  is  no  single-valued 
function  y  =  y(x)  describing  C  in  a  neighborhood  of  x  =  x, . 
The  curve  C  must  therefore  be  concave  toward  the  fluid  in  a 
neighborhood  of  P,,  i.e.,  locally  the  curve  will  lie  to  the 
left  of  X  =  x^  (see  Fig.  6). 

Suppose  we  move  a  vertical  line  toward  the  curve  C,  from 
the  left.   Then  this  line  must  be  tangent  to  C  at  a  point 
Pg^-^p^y?^  where  x^  *=  x,  .   Note  that  if  s,  and  s^  are  respec- 
tively the  values  of  s  corresponding  to  P,  and  P^,  then  we  may 
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take  Sp  >   s,.   Near  tl        Pp,  however,  we  see  that  the 
angle  th      ent  makes  with  the  positive  x-axls  Is  decreasing 
to  7r/2.   This  means  i^  is  negative  In  a  neighborhood  of  Pp. 
This  contradicts  the  Inequality 

at  Pp,  which  says  that  k   must  be  positive.   Hence  the  curve  C 
can  have  no  vertical  tangent  at  any  point  and  has  a  non- 
parametric  representation  y  =  y(x),  where  y(x)  is  a  single- 
valued  function  of  x. 
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CHAPTER  4 
THE  SECOND  VARIATION  OF  M 

In  this  chapter  we  set  up  a  variational  principle  for  the 
wave.   This  principle  asserts  that  If  A  and  E  are  held  fixed, 
then  the  stream  function  ^  describing  the  wave  minimizes  M.   We 
recall  that  In  Chapter  3  we  showed  how  we  could  restrict  the 
class  of  competing  functions  which  would  enter  Into  an  existence 
proof  for  waves  by  variational  techniques.   If  we  limited  our- 
selves to  one  period  of  the  free  surface,  then  symmetrlzatlon 
indicated  that  the  wave  had  the  non-parametric  representation 
X  =  x.(y)  with  just  two  branches,  which  are  symmetric  about  a 
trough  or  crest,  and  Hopf ' s  lemma  implied  that  the  wave  had  the 
non-parametric  representation  y  =  y(x).   Thus  the  class  of 
admissible  functions  could  be  limited  to  functions  having  these 
two  properties. 

We  note  that  only  certain  values  of  E  will  give  a  solution. 
This  is  because,  if  we  fix  A  alone,  then  a  minimum  is  furnished 
for  E  when  we  consider  the  flat  case  of  uniform  flow.   Thus  E 
must  be  larger  than  this  value.   If  we  allow  the  amplitude  to 
get  too  large,  then  again  we  presumably  get  no  solution  because 
the  wave  must  break. 

We  show  in  this  chapter  that  the  above  variational 
principle  is  correct  by  deriving  a  formula  for  the  second  varia- 
tion of  M  and  showing  it  to  be  strictly  positive  when  A  and  E 
are  held  fixed  up  to  second  order. 
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We  have  already  shown  in  Chapter  2  that  if     .  crlbes 
a  wave  then 


If  also 


and 


then  clearly 


6F  =  5(M  +  (A-c^)A  -  E)  =  0  . 


5A  =  0 


6E  =  0 


6M  =  0  . 


2 

Let  us  denote  the  second  variation  by  6  .   We  will  establish 

the  following  result. 
Theorem.   If 


and 


BA  =  B^A  =  0 


5E  =  B^E  =  0 


then  not  only  5M  =  0,  but  also 


5^M  >  0 


provided  tJ/   describes  a  wave.   This  will  imply  that  the  minimum 
problem  is  correct  for  waves.   In  this  paper,  however,  we  are 
able  only  to  prove  the  theorem  for  waves  of  sufficiently  small 
amplitude. 


k2 


4. I. a.  Derivation  of  the  Second  Variation  of  M 


Let 


J(e)  =  j(/  F(x,y,^,^^,ty)dA 


D 

be  a  functional  which  varies  analytically  with  the  parameter  e, 
By  Taylor's  theorem,  we  may  write  the  series 


J(e)  =  J(0)  +  (g)e  +  I  (^)e^+  ...  . 

de 


Then  the  first  variation  of  J  is  defined  as 


6J  =  (§|) 


£=0 


and  the  second  variation  is  defined  as  the  next  term  in  the 
series,  that  is. 


p2^   1  ,d^j, 
^  ^   =  2    ^dF-) 


e=0 


In  Section  2. II. a.,  we  derived  the  formula 


6M 


-   -        (c  -  q  )5nds 


for  the  first  variation  of  M.  We  took  n  to  be  the  inner  normal 
to  the  analytic  curve  C  which  forms  the  upper  boundary  of  D  and 
we  let 

6n  =  ep( s) 


^3 


whe;      )  Is  a  real  analytic  function  of  the  arc  length 
parameter  s  ^^-..t;  C. 

In  this  section  we  derive  the  following  formula  for  ST*!: 

S^M  =  -  ^J   (c^  +q^)K(5n)^ds  +  JJ   (V(6V))^dA  , 


where  n,  C,  and  D  have  the  same  meanings  as  described  above  in 
the  formulation  of  the  expression  for  5M,  and  k  is  the  curvature 
along  C.   We  let  tj/  be  a.   stream  function  defined  on  D.   As  before 

V  =  0 

on  C, 

V(0,y)  =  V'd^y) 

for  all  y,  and 

{V^)2  =  0(1) 

as  y  -»■  -00  . 

We  obtain  the  domain  D*  by  varying  only  the  upper  boundary 

of  D.   Let  C*  be  the  upper  boundary  of  D*  and  let  f*   be  the 

corresponding  stream  function  defined  on  D  ,  which  has  the 

following  properties: 


throughout  D*, 


for  all  y,  and 


as  y  -»■  -00  . 


^Tf/*     =     0 

V*(l,y)  =  V'*(0,y) 
(VV'*)^  =  0(1) 
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We  make  the  assumption  that  C  and  C*   have  the  same 
Initial  and  terminal  points.   Without  loss  of  generality  the 
formula  for  the  second  variation  of  M  will  then  be  derived 
under  the  additional  hypothesis  that  C  and  C*  intersect  each 
other  once  between  their  common  end  points  (see  Pig.  7).   The 
techniques  Involved  in  proving  the  formula  for  the  case  of  one 
intersection  are  easily  extended  to  an  arbitrary  number  of 
intersections. 

Let  C  =   C-^  +  C^  and  C*  =  C*+C2,  where  C,  and  C*  bound 
D-D*  and  Cg  and  C*   bound  D*-D.   Since  C  is  an  analytic  curve  and 
p(s)  is  a  real  analytic  function,  C  is  also  analytic.   This 
implies,  by  an  application  of  the  Schwarz  reflection  principle, 
that  ^  and  f*   can  be  extended  outside  of  their  original  domains 
of  definition  in  such  a  way  so  that  they  remain  harmonic.   Hence 
we  can  take  C  and  C*  so  close  that  ^  and  ip*   are  harmonic  in 
D-D*  and  D*-D. 

Let 

5^  =  f*  -  Tp    . 

Then  bip   is  harmonic  in  (D-D*)  U  (D*-D)  U  (Dn  D*),  and  the  normal 
derivative  of  6^  vanishes  on  the  vertical  boundaries  of  D  and  D*. 
We  will  need  the  expansion  of  6t//  for  y  large  and  negative.   From 
Theorem  2. I. a.,  we  have 

00      p 

■d/  =  y    +  cy  +  S        a  e  '^^^   cos  ( 2Tmx  +  h) 
^        ■'  o        "^        ^ — =1-  n  n 

n=l 

and 
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i>*  =  y*+cy  +  y~  a*e^"y  cos   (27rnx  +  ^*) 
°  n=l 

near  y  =   -oo  .      Thus 

00  Q  ~ 

5^  =  y!  -y^  +  ZZ  ^n^  ^°^   ^^'^  +^n^ 

n=l 

From  this  expansion,  we  see  that 


<   Ke^y 


If  y  is  large  and  negative.   The  expansion  for  6^  also  implies 
that  the  Dirichlet  integrals  for  5^  over  D  and  D  both  converge. 
In  Section  2.II.b.,  we  saw  that  we  may  write 

(M*-M)  -  c^(A*-A)  =  JJ    V^-  V^*dA  -  JJ    V^.  V^*dA  . 

D-D*  D*-D 


Since  V*  =  ^  +  5^, 


//     VV^-  V^*dA   -      //     VV/ .  V^*dA   =      //     V^.  V(V/+6V/)dA 


D-D*  D*-D  D-D* 


-  //   V^/  .  V(V  +  5V')dA 


D*-D 


2, 


Substituting  (M*-M)  -  c^(A*-A)  for  the  right  hand  side  of  the 
equation  and  expanding  the  left  hand  side  gives 
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(M*-M)  -  c^(A*-A)  =   rr  (V^)^dA  +   /yV^-  V(6t)dA 


D-D*  D-D* 


D-D*  \  D*-D 


+   //  V^.  V5T//dA  +   //  (V(5i//))^dA 
D*-D  D*-D 


An  application  of  Green's  theorem  to  the  domains  D-D*  and  D*-D, 
using  inward  normals  along  the  boundary  curves  with  respect  to 
the  domains  D  and  D*,  respectively,  gives 

(M*-M)  -c^(A*-A)  =  JJ  {\/f)^dk  +J    f  -^   (5^)ds 

D-D*  C* 

-  f  "P  ^   (6^)ds  -  f  Jf  {yip)^dk 

C,  \  d"*-d 


-/*^-^/*^- 


C|  Cg 


Since  f   is  zero  along  C  and  hence  along  C^  and  Cp,  and  since 
Tp   =  -6^  on  C*,  by  combining  the  two  non-zero  line  integrals,  we 
have 

(M*-M)  -  c^(A*-A)  =  JJ  {Vf)^dA    -     ^  (V^)^dA  -^(5^)  ^  (5t)ds  . 

D-D*  D*-D  C* 

Let  Df  be  the  domain  which  has  as  its  upper  boundary  the 
curve  C*,  as  its  lower  boundary  the  line  y  =  -k  where  k  >•  0, 
and  whose  side  boundaries  are  the  usual  vertical  lines.   Then 
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An  application  of  Green's  theorem  to  the  Inte- 


r<.  ^^    V*. 


gral  -J  J^  (5^)ds 


C* 

0 
-  f  i^i')   -^   (5^)ds  =^V(5^)^dA  +  f     (5^)  ^  (5V)ds 

k  x=0 

-k  0 

+  /  (5^)  -^   (5^)ds  +  J     (5^)  ^  (5V')d£ 

0  1 

x=l  y=-k 

for  every  k  >  0. 

Since  we  know  that  5^  is  periodic,  the  second  and  third 
Integrals  on  the  right  hand  side  of  the  equation  cancel  each 
other.   It  was  shown  above  that  -r—  (5^)  -♦  0  exponentially  as 
k  -►  OD  .   This  implies 


11m   /   (5V)  -^   (5V')ds  =  0  . 
k  —  OD  -^       °" 


0 
y=-k 


Thus  we  have,    by  taking  the  limit  as  k  -►  co  , 
-/  (BV')  ^   (5V)cls  =  ff  [V{bi,)faA 


C*  D* 


(M*-M)  -  c^(A*-A)    =     ^  (V^)^dA    -      |J'(VV/)^dA    +      /T  [  V(  bV')  ]^dA 

D-D*  D»-D  D* 

Replacing  the  Integral  over  D*  with  an  integral  over  D  gives 
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(M*-M)  -  c^(A*-A)  =  J^(V^)^dA  -   /T  (V^)^dA 


D-D*  D*-D 


+  ff   [V(5^)]^dA  +  O(e^) 


D 


We  now  express  dA  In  terms  of  normal  and  tangential 

2      2 

coordinates  and  replace  (V^)   by  q  to  get 


(„..„,  .c2(A*-A)  =       ff  q^aU^,^,^ 


(D-D*)-(D*-D) 

+  JJ   [V(5^)]^dA  +  0(£^)  . 


D 


Inserting  an  expression  for  the  Jacobian  valid  up  to  second 

2 
order  terms  and  writing  q   in  terms  of  its  value  on  the  curve  C 

through  Taylor's  theorem,  we  get 

(M*-M)  -c^(A*-A)  =       JJ  (q^  +  |^nj(l+Kn)dnds 

(D-D*)-(D*-D)^        / 

+  jf   (V(5^))^dA  +  O(e^)  . 

d" 


2 
Here  K  denotes  the  curvature  on  the  carve  C  and  q  denotes  the 

2 

value  of  q  along  C.   Integrating  with  respect  to  the  normal 

coordinate  and  dropping  off  terms  of  higher  than  second  order 
gives 
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(„*_„).  c2(A»-A)  =/q2(6n)ds+/  |^  ^  d 


2  ,..x2 

"3 


+ 

c 


f^   q^(6n)^ds  +ff  [V(6V)]^dA  +  O(e^) 


Since  up  to  second  order 

A*-A  =  -  J  [1+^   6n)(6n)ds 


-x_2        p 

and  since  ■?3_  =  _  2Kq  ,  we  get  the  following  result; 
an 

M*-M  =  -  f  {c^  -q^)bnds   -  i  j  ( c^  +  q^)K(  &n)^ds 


+ 
D 


fj   [V{6^)]^dA  +  0(6^) 


Because  the  second  variation  of  M  consists  of  all  the  above 
terms  that  are  of  second  order  in  e,  we  have 

S^M  =  -  ^J   (c^+q^)K(5n)^ds  +  JJ   [V(6V')]^dA  . 
C  D 

U.I.b.  The  Second  Variation  of  M  For  a  Wave 

In  Chapter  2,  we  showed  that  if  the  area  A  and  the  poten- 
tial energy  E  were  held  fixed,  and  if  f   described  a  wave,  then 
the  first  variation  of  M  vanished.  We  now  investigate  the 
second  variation  of  M, 
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5^M  =  -  IJ  (c^+  q^)l<(5n)^ds  +  JJ   [V(5^)]^dA  , 


D 


when  these  conditions  are  imposed. 

If  A  and  E  are  held  fixed,  then  we  may  say  that 


and 


5A  +  6^A  =  0 


&E  +  5^E  =  0 


up  to  second  order  in  e.   It  is  easily  seen  that  we  may  write 
5A  +6^A  =  -^  (1+1  &n)6nds 


and 

2gy  5nds  -J    (gyK  +  g  ^)  (  5n)^ds 


5E  +5^E  =  -  /  2Ky  5nds  -  /  f  evK  +  c-  ^U  ftn  l^r 


Consider  the  following  sum 

5^M  =  (5M  +  5^M)  +  (A-c^){6A  +5^A)  -  (5E+5^E) 

-J    (-A  +q2  +  2gy)6nds    +  [^    ( -A  +  q^  +  2gy  )k(  6n)^ds 
C  ^C 

+  /^ll^  +  g|^)(6n)2ds    +j^[V(6V/)]2dA    +  0{e^)    . 
C  D 

Since  ^  is  a  wave,  it  satisfies  Bernoulli's  law  on  C,  i.e., 

q^  +  2gy  =  A  . 
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This  shows  that  the  first  and  second  integrals  In  the  above 
equation  vanish.   Thus  we  arrive  at  the  formula 

6^M  =  -^/^  (q^  +  2gy)(6n)^ds  +  ff  [V(5V/)]^dA  . 


By  Hopf ' s  lemma  for  subharmonic  functions 

^  (q^  +  2gy)  <  0  . 

2 
This  implies  that  the  first  integral  in  the  formula  for  6  M 

is  always  negative.   Note  that  the  second  Integral  is  always 

positive. 


4. II.   Investigation  of  the  Sign  of  the  Second  Variation  of  M 

We  now  have  an  expression  for  bT^i  in  the  case  of  a  wave 
when  variations  are  made  holding  A  and  E  fixed.   This  expression 
is 

6^1  =  -^/^  (q^+2gy)(6n)^ds  +  ff  [^{bf)fdA    , 
C  D 

where  the  first  term  is  negative  and  the  second  term  is  positive. 
We  want  to  show 

6^  >  0 

subject  to  two  linear  constraints  on  6n, 

5A  =  0  ,    6E  =  0  . 
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We  note  that  showing  the  above  inequality  is  equivalent  to 
showing  that  the  quotient 


ff  [V{b^|^)fdA 


^  >  1 


f^   (q2+2gy)(6n)2ds 


1  r  S  ,_2 

C 


subject  to  the  constraints 

5A  =  -  /  5nds  =  0 
and 

5E  =  -  /  y  5nds  =  0  . 

Once  we  have  proved  this,  we  may  conclude  that  the  wave  mini- 
mizes M  for  fixed  A  and  E  and  the  variational  principle  we 

stated  for  the  wave  is  correct. 

2 

First  we  will  consider  the  corresponding  expression  5  M 

for  the  flat  case  of  uniform  flow  and  show  it  to  be  strictly 
positive.   Then  we  will  deduce  from  this  result  that  dtA   is 
strictly  positive  for  the  wave  provided  we  assume  that  we  are 
close  to  the  flat  case. 

4. II. a.  The  Sign  of  the  Second  Variation  of  M  for  the  Flat  Case 
Recall  that  for  a  wave  we  have 

5^  =  -1/  (|^  +  2g  |^){6n)2ds  +ff   [V(6^)]^dA 
C  D 
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when  A  and  E  are  stream  function  f      flat 

case  is  glvtw  . ,. 


T^^MS,  we  may  rewrite  the  above  expression  for  ST'I  in  the  form 


0 


where  R  is  the  semi-infinite  rectangle  0  <  x  ;^  1,  y  j^  0. 
Through  an  application  of  Taylor's  theorem  and  noting  that 
q  =  |°3l|  on  C,  we  may  write 


q  c 

LStituting  this  result  for  (6n)  ,  we  see  that  showing 
:■•  >  .  is  equivalent  to  proving  that 

-^. >  1 


0   ° 


(6l^)^dx 


the  two  constraints  that  A  and  E  are  held  fixed. 

o 
In  order  to  evaluate  g/c  and  to  investigate  our  con- 
straints, we  consider  the  flat  case  as  the  limit  of  the  classi- 
cal linearized  soluti      whlnh  the  amplitude  h  is  small.   We 
the..,.....,  .,^.,.^3  this  1.,.-.-.  ...d  case  briefly  (cf.  \9]). 


5^4 


Let  ^(x,y)  be  the  velocity  potential.   For  the  linearized 
case  we  may  redefine  ^  by  the  expression 

J(x,y)  =  ex  +(t)(x,y)  , 

where  the  gradient  of  ^   will  be  inf initesimally  small.   We  may 
write 

Vf  =    (c,0)  +(^^,<t)y)  , 

where  (({)  ^if)  )  exhibits  the  velocity  due  to  the  wave  motion. 
X   y 

We  now  let  y  =  r[{x)    be  the  equation  of  the  free  surface 
of  the  wave,  where  r)(x)  will  be  as  small  as  we  wish.   We 
substitute  •n(x)  and  (V^)   into  Bernoulli's  equation  to  get 

H  =  gTi  +  2 —  > 

where  H  is  a  constant.   Expanding,  we  obtain 

c^+  2c(t)  +  <j)^  +  (f)^ 
H  =  gri  +  1 ^ ^  , 

where  the  terms  on  the  right  hand  side  of  the  equation  are 
evaluated  at  y  =  r)(x). 

Since  t]  and  grad  ^   are  small,  we  may  linearize  the  above 
by  simply  dropping  the  quadratic  terms  and  expanding  in  a 
Taylor's  series  about  y  =  0.   Thus,  applying  Taylor's  theorem 
to  Bernoulli's  equation,  we  have 


IT    2    ci 
He     ^x 


^  =  77  -  077 


g    2g    g 
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y=o 


Now  we  establish  the  klnetl  r  the  free  surface 

_, ^  ...,   ime  derivative  of  the  equation  y-'n(x)  =  0  and 

set'     t  equal  to  zero,  which  gives 

ly  -  ^xlx  =  °  • 
Linearizing  the  equation  gives 

or 

'x    c 

From  the  equation  for  t^(x)  which  was  derived  from  Bernoulli's 
law,  we  have 

'x     g  ^xx 
Combining  then  two  equations  for  r\^,    we  get 

'x     g  ^xx    c 


which  implies  the  classical  result 

2 
^y     g  ^xx 


The  linearized  equation  Just  derived  has  the  family  of 

solutions 

(|)  =  A  e*^  sin  kx  , 
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where  k  =  g/c  .   Our  periodic  boundary  restriction  gives  us 


k  =  2Tn,    and  hence 


27rn  =   -^  . 


2 

We  shall  use  this  relationship  in  the  expression  for  5  M  in  the 

flat  case. 

We  combine  the  equations  for  ^   and  make  the  following 
substitution  to  display  the  small  magnitude  of  the  wave: 

2 

h  =  . 

g 

We  have 

T]   -   h[l-cos    27mx ] 
and 

(t)   =   che        '^cos   2Tmx    . 

We  further  restrict  ourselves  by  taking  the  solution  which 
has  period  one.   Then 


ri  =  h[  1  -  cos  2Trx]  , 

^  =  cy  +ch  e  ^^  cos  27rx  , 

and 

-%  =  27r  . 
c 

2 

We  may  now  write  down  the  expression  for  5  M  for  the  flat 

case  by  considering  it  as  the  limit  as  h— •■  0  of  the  linearized 
solution  of  period  one.   This  is 
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1 

r 

0 


1 
e^M  =  -27r  r  (6V')^dx  ^  ''  [V(6V)]^dA  , 


provided  A  and  E  are  held  fixed.   Thus  we  wish  to  show  that 

^[V(6^)]2dA 


^-. >  1 


2Tr  /  ibf) ^dx 


0 

subject  to  two  linear  constraints  on  bTp. 

We  must  investigate  the  variational  constraints.   As 
mentioned  before,  we  assume 


.1 

A*-A  =  /   6ndx  =  0  , 
0 


-=/ 


where 


Thus 


6n  =  -  -^  (6V/)  . 


1 

/   6i^dx  =  0  . 
0 

Similarly 

5E  =  -  /  y  6nds  =  0  . 
C 

To  derive  from  this  a  nontrivial  constraint  in  the  case  of 
uniform  flow,  we  consider  the  limit  of  the  linearized  wave  with 
free  boundary 

r]   =  h(l  -  cos  27rx)  . 
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We  may  write 

1 


/   h(l-  cos  2Trx)B^dx  =  0 


0 

for  each  h.   Using  the  earlier  constraint  on  A,  we  conclude 

that 

1 
/   (cos  27rx)5^dx  =  0  . 
0 

Our  problem  "becomes   to  minimize   the   quotient 
I   =-^. 


2ir\i  dx 


0 

y=o 


over  the  class  S  of  functions  u  defined  over  the  semi-infinite 
rectangle  R  such  that 

Au  =  0 

throughout  R, 

u(0,y)  -  u(l,y) 

for  all  y  _^  0,  and  u  -»■  0  exponentially  as  y  -^  -oo  .   As  a  result 
we  are  led  to  the  eigenvalue  condition 


^  +   27rAu  =  0 
dn 
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on  the  Interval  (0.1)  along  the  x-axls.   The  corresponding 

F-  •  v-T  :  .  c  igenfu!.-- '.  -  ..s  are 


1 

u_  = 


and 

1   kv 
u,  =  —  e   cos  kx  , 

^  A 

where  k  =  27r^  and  X  =  1,2,.-.  •   Thus  zero  is  the  lowest  Stekloff 
eigenvalue  and  the  remaining  X's  are  the  other  Stekloff  eigen- 
values.  These  are  related  to  the  linearized  waves  discussed 
earlier. 

The  elgenfunction  u^  gives  the  minimum  of 


1=  ^ 


If  ^^^^^^^ 


1 

^ 

0 


/ 


2Tru  dx 


over  all  functions  in  the  class  S.  This  minimum  value  is  of 
course  zero.  The  elgenfunction  u,  minimizes  the  above  ratio 
over  all  functions  in  S  orthogonal  to  u  in  the  sense 


1 

udx  =  0  . 
0 


/ 


We  note  that  6V'  satisfies  this  orthogonality  requirement  because 

1 

/  bi/dx  =   0  . 
0 
We  have 
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J]  <%'"" 
^. =1 , 


and   thus 


27ru,  dx 


ff  [V(5^)]^dA 

■4 >    1 


27r(5^)^dx 


0 


To  get  the  next  eigenvalue,  we  minimize  the  above  ratio 

with  respect  to  all  functions  in  S  orthogonal  to  both  u  and  u, . 
^  *="  o      1 

A  function  in  this  class  has  the  following  properties: 

1 
/   udx  =  0 

0 

and 

1 

/   u  cos  27rx  dx  =  0  . 

0 

Since  the  next  eigenvalue  is  A  =  2,  for  all  such  functions  we 
have 


[f   (Vu)2dA 


R     >  2 


1 


27ru  dx 


Because 


0 

i 

/   &^  cos  27rxdx  =  0  , 
0 
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that 


JJ  [y{b^)fdk 


-5-^j >  2  >  1  , 

'   2Tr  (  5V )  ^dx 
0 


Is  what  we  set  out  to  prove. 

Our  analysis  shows  that  in  the  flat  case  of  uniform  flow 
we  have 

6^  >  0 

under  the  single  constraint  5A  =  0;  but  when  both  constraints 
6A  =  5E  =  0  are  imposed  we  obtain  the  stronger  inequality 


)^ 


>  0 


The  latter  condition  will  enable  us  in  the  next  section  to  esti- 
mate 5'T*'!  for  waves  of  small  finite  amplitude. 

.  ;  .  The  Sign  of  the  Second  Variation  of  M  for  Small  Amplitude 
Waves 

Let  us  assume  that  our  wave  is  close  to  the  solution  for 
the  flat  case  described  in  the  previous  section  (4. II, a.)-   The 
stream  function  V  niay  be  expanded  in  a  power  series  In  the  ampli- 
tude h  which  approaches  the  flat  solution  as  h  -►  0.   We  have, 
therefore, 

uy  +  h^^  +  h^Vg  +  •  •  • 
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for  h  sufficiently  small.   In  this  situation  we  will  show  that 

B^  >   0 

provided  A  and  E  are  held  fixed.   The  result  Implies  that  the 
wave  minimizes  M  for  fixed  A  and  E. 
We  have  the  expression 

e^M  =  -  r  p(B^)^ds  +  I    [V(5^)]^dA  , 

C  D 

where 

q 

and 


2 

We  note  that  p  >■  0.   The  expression  for  5  M  will  be  positive 

provided  we  can  show  that 


J]    [V(6^)]2dA 

^ >   1 


2 
p(5^)  ds 


when  two  linear  constraints  arising  from  fixing  A  and  E  are 
imposed  on  6^^. 

Let  us  consider  these  constraints.   The  constraint  on   A 
implies 

5A  =  -  /  Bnds  =   -  -  bi/ds   -  0  , 
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■Ives  us 


/aaa=  /M.x.o(h: 


or 


C  0 

y=0 


J. 
I      bTl/dx  =  0(h)    . 


0 
y=o 


The  constraint  on  E  can  be  written  as 

5E  =  -  /  y  5nds  =  0  . 
C 

Combining  these   constraints  gives 


J. 
0  =  ^   /  y   5nds  =  ^  (cos   2Tnc)6tdx  +  0(h)    , 


0 
y=0 


which  implies 


X 

/   (cos  27rx)B^dx  =  0(h)  . 


0 

y=o 


Let      troduce  the  Nuemann's  function  so  that  we  may 
replace  ^.wc  .   '  '  ent 

ff   (Vu)^dA 
I  =  ^ 


fpu 


^ds 


by  an  equivalent  quotient  of  two  quadratic  forms  defined  over 


6k 


the  same  space.   The  Neumann's  function  N(z,C)  with  respect  to 
a  domain  D  is  a  normalized  fundamental  solution  of  Laplace's 
equation.   There  are  several  ways  to  introduce  the  Neumann's 
function.   For  the  application  that  we  have  in  mind  it  Is 
convenient  to  define  the  normal  derivative  in  the  following  way: 
On  C,  we  take 

SN  ,  ).>,        2ira. 
(z,C)  = 


^  v->s^  -  -^  > 


where 


q      a  ' 


C 


and  on  the  vertical  boundaries  of  D 

an 

This  choice  of  normal  derivative  is  compatible  with  the  condition 


1     r    SN  ,  ... 


which  holds  because  of  the  pole.   For  harmonic  functions  u  which 
satisfy  the  condition 

c 

this  choice  of  --r—  simplifies  the  representation  we  derive  for  u 
in  terms  of  its  normal  derivative  on  SD.   The  condition  on  the 
area 
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-  ds   =  0 

■1 
C 


Is   such  a  function.      In  particular,    we  may  write 
ippllcation  of  Green's   theorem 

u(a    =  -^  f  [u{z)  ^  (z,0    -   N(z,C)  ^   (z)]ds 


^/      N(z,C)^a3./^d3. 


This  reduces  to 


u(u  =  -^/N(z,a|^ 


ds 


provided  we  restrict  ourselves  to  the  class  S  of  harmonic  func- 
tions u  which  have  zero  normal  derivative  along  the  vertical 
boundaries  of  S  and  such  that 


/! 


■^  ds  =  0  . 
C 

We  will  also  assume  that  functions  in  s'  decay  exponentially  to 
zero  as  y  —  -oj  . 

Using  the  Neumann  function  N(C,z)  for  the  domain  D,  we  may 
rewrite  the  quotient 


1 


rr  (vu)^dA     y^ui^ds 


/  pu  ds         /  pu  d 


s 
C 
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by  applying  Green's  theorem  and  using  the  above  representation 
for  functions  u  in  S'.   We  then  get 


27r//N(C,z)  1^  (z)  |H  (Odsdt 


^  C  C 


y'/|pN(C,t)N(t,z)  |H  (z)  |H  (Odsdtda 

c  c  'c 

We  introduce  a  Hilbert  space  S*  which  Is  a  subspace  of  the  spj 
Lp  of  Lebesgue  square  integrable  functions.   The  subspace  is 
defined  by  the  orthogonality  condition 

J  J    M^  w(z)dsdt  =  0 

c  c 

for  all  w(z)  in  S  .   This  condition  is  equivalent  to 


Mds  =  0 


when 


Let 


q 

C 


w(z)  =  -^  (5^) 


K(C,z)  -  I    pN(C,t)N(t,z)dt 
C 


Then  minimizing  I  over  functions  in  S   is  equivalent  to  mini- 
mizing the  quotient  Q.,/Qp  of  the  two  quadratic  forms 


Q3_  =  J  J  N(C,z)w(z)w(adsdt 


C  C 

and 
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i.,       K(C,z)w(z)w(C)dsdt 

^       _  ■_' 

C  C 

over  the  functions  In  S*. 

2 

We  now  proceed  to  show  the  desired  Inequality  6  M  >■  0  for 

fixed  A  and  E  by  comparing 

/  /  N(C,z)w(C)w(z)dsdt 
I  =  ^  ^ 


y^J  K(C,z)w(z)w(C)dsdt 
C  C 

with  the  corresponding  expression  in  the  flat  case  of  uniform 
flew.  Specifically  we  show  that  for  h  sufficiently  small  and 
0  <  6  <  1 


r  rN(C,z)w(z)w(C)dsdt 


_2_2 >  2  -  6  >   1 

r  rK(C,t)w(z)w{adsdt 
C  C 

for  functions  w  In  S*  which  satisfy  the  following  condition 
arising  from  fixing  A  and  E: 

J'y  K(^,z)w(z)u^(C)dsdt  =  0(h)  , 
C  C 

where  u-,  =  e*^^  cos  Sttx  is  the  first  Stekloff  eigenfunction 
discussed  in  'I. II. a.   This  Inequality  involving  the  Neumann's 
function  will  give  us  ' 


f.9. 


6^M  >   0 

for  fixed  A  and  E. 

For  the  flat  case  it  is  easily  seen  that  for  functions  in 
s',  the  quotient  I  may  be  rewritten  in  the  form 

1   1 

//  No(-'^)|S(^)|S(-)^^^- 

^     111 

f   f   f     2Tr[N^(t,x)N^(t,C)  1^  (C)  1^  (x)]dtdsdx 
0   0   0 

where  N  (x^C)  is  the  Neumann's  function  for  the  sem.i-inf inlte 
rectangle  R.   Then  except  for  the  eigenfunction  u^  =  l//27r,  the 
Stekloff  eigenf unctions  described  in  Section  4. II. a.,  when 
evaluated  at  y  =  0,  are  also  the  eigenfunctions  for  the  above 
quotient  I  of  quadratic  forms.   The  eigenvalues  are  unchanged. 
Thus  we  have 

1   1 


J   J    N^(x,C)w{x)w(C)dxds 
f    f    K^(x,C)w(x)w((;)dxds 


0  0 
where  w(x)  is  in  S*  and 


K^(x,C)  =  j   27rN^(x,t)N^{t,s)dt 
0 


The  orthogonality  condition 
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1 

COS  27rx  dx  =  0 
0 


/" 


translates  through  the  equality 

1  1 


/     u  cos  2Trx  dx=/      /     K(x,t)4^(t)    cos   27rx  dxdt 


0      0 


into   the   constraint 

1      1 
/      /      K   {x,t)w(t)    cos   27rx  dxdt  =   0    . 
0     0 

Imposing  this  contralnt  will  therefore  [cf.  4. II. a.]  give  us 

1  1 
J  J    NQ(x,aw(x)w(C)dxds 

1  1 -  ^  • 

'^     ^    K^(x,aw(x)w(C)dxds 

u  0 

We  note  that  x.ne   constraints  for  the  non -linear  case  differ 
from  the  constraints  Imposed  in  the  flat  case  of  uniform  flow 
by  0(h).   Since  p  =  27r  +  0(h)  and  since  the  Neumann's  function, 
an^  '    -  '-'- '   ^--.---.    -  r.   ^j^g  ^^Q  quadratic  forms  Q,,  Qp  change 
continuously  with  the  domain,  we  see  that  for  h  sufficiently 
small 


/  /  N(C,z)w(z)w(/;)d8dt 


^  ^  >2-6>l, 


j  y  K(C,z)w(z)w(adsdt 
C  C 
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where  5  *=  1 .   Since  -r—  (5^)  is  in  S  and  satisfies  the  constraint 

y^y'K(z,Ow(z)u^(C)dsdt  =  0(h)  , 
C  C 

we  may  assert  that  the  above  inequality  holds  for  w  =  -^   (6?//). 
Hence  we  conclude  that 

6^  >   0 

for  a  wave  provided  the  amplitude  is  sufficiently  small,  and  L 
and  E  are  held  fixed. 
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